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Abstract We consider the solution u : [0,°°) x % d — s- [0,°°) to the parabolic Ander- 
son model, where the potential is given by (t,x) i— > ydy, (x) with Y a simple sym- 
metric random walk on If 1 . Depending on the parameter 7 e [— °°,°°), the potential 
is interpreted as a randomly moving catalyst or trap. 

In the trap case, i.e., 7 < 0, we look at the annealed time asymptotics in terms of the 
first moment of u. Given a localized initial condition, we derive the asymptotic rate 
of decay to zero in dimensions 1 and 2 up to equivalence and characterize the limit 
in dimensions 3 and higher in terms of the Green's function of a random walk. For a 
homogeneous initial condition we give a characterisation of the limit in dimension 1 
and show that the moments remain constant for all time in dimensions 2 and higher. 
In the case of a moving catalyst (7 > 0), we consider the solution u from the per- 
spective of the catalyst, i.e., the expression u(t,Y t +x). Focusing on the cases where 
moments grow exponentially fast (that is, 7 sufficiently large), we describe the mo- 
ment asymptotics of the expression above up to equivalence. Here, it is crucial to 
prove the existence of a principal eigenfunction of the corresponding Hamilton oper- 
ator. While this is well-established for the first moment, we have found an extension 
to higher moments. 
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1 Introduction 

The parabolic Anderson model (PAM) is the heat equation on the lattice with a 
random potential, given by 

f t u(t,x) = KAu(t,x) + %(t,x)u(t,x), (t,x) G (0,°o) x Z d , 

u(Q,x) = uq(x), x <E Z d , 

where K > denotes a diffusion constant, uq a nonnegative function and A the 
discrete Laplacian, defined by 

A f (*):= E \f(y)-f( x )]> xeZ d ,f:Z d ^M. 

y ez d : 

|.v-v| = l 

Furthermore, t, : [0,°°) x If 1 — > K is a space and time dependent random potential. 
We deal with the special case that the potential is given by 

S(t,x) = ySr,(x), (t,x)e {0~)xZ d , 

with a simple symmetric random walk Y with generator pA that starts in the origin 
and a parameter 7 G [— °°, °°) called coupling constant. In this paper we analyse the 
large time asymptotics after averaging over the potential which is usually referred 
to as annealed asymptotics. We denote expectation with respect to the potential £ 
by (•). 

One possible interpretation of this system arises from chemistry. Here, u(t,x) de- 
scribes the concentration of reactant particles in a point x at time t in presence of 
a randomly moving particle. In the case 7 < 0, the particle acts as a decatalyst (or 
trap) that kills reactant particles with rate — 7 at its position. In the case of positive 7, 
we consider a catalyst particle that causes reactants to multiply with rate 7. In both 
cases (u(t,x)) is interpreted as the averaged concentration. For further interpreta- 
tions and an overview over the PAM see for instance HGM90I . 1CM94I . IIM94I and 
IGK051 . 

Annealed asymptotics in the case of a positive coupling constant 7 have already 
been investigated in 1GH061 . In the present work, we derive similar results with re- 
gard to the expression u (t,x) := u (t,Y t +x), which can be interpreted as the particle 
concentration in a neighbourhood of the catalyst. In addition to logarithmic asymp- 
totics in terms of Lyapunov exponents, we derive asymptotics up to equivalence for 
most of the parameter choices where exponential growth is observed. 
The case that 7 is negative has to the best of our knowledge not been investigated so 
far. Its analysis relies on techniques quite different from those in the catalyst case as 
a functional analytic approach proves unfeasible here. We calculate moment limits 
dependent on the model parameters and, in the case of moment convergence towards 
zero, specify the convergence speed up to equivalence. 

Whereas the PAM with time independent potential or white-noise potential is well 
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understood, some other time dependent potentials have just been examined recently. 
In MGdH06L BGdHMl IB and I1KS03I , for instance, the authors investigate the case 
of infinitely many randomly moving catalysts. In IICGM1 111 the authors deal with 
the case of finitely many catalysts, whereas the article j DGRS 1 1 j] is dedicated to 
a model similar to the case of infinitely many moving traps. Further examples of 
time dependent potentials can be found in IIGdHM071 , HGdHM09al , IGdHMlOL 
IIMMS11I . and the recent survey IIGdHM09bll . Within these proceedings, HKSllL 
ILMl 111 and BMZl 111 deal with the parabolic Anderson model with time-independent 
potential. 

In Section 12.11 we analyze the PAM with localized initial condition mo = 8 Z and 
y < 0. Let 

M z (t):=( £ u(t,x)), (t,z)€[0,<x>)xZ d , 



denote the expected total mass of the system at time t if the solution is initially 
localized in z and the trap starts in the origin. We find 

Theorem 1. For d = 1 , 2 and every z £ Z° ', 

(0 M z (t)~^^ri, t^oo far d = l; 

(ii) M z (t) -47r^(logf) _1 , f^°° for d = 2, 

and 

Theorem 2. For d>3 and every z £ Z d , 

limM-(f) = H - r^Gi(z), 

where G K denotes the Green's function of a random walk with generator kA. 

Remark 1. Theorems Q] and |2] can be generalized to all initial conditions with com- 
pact support without much effort. 



In Section I2T21 we analyze the case of a homogeneous initial condition uq = 1. We 
find that in dimensions 2 and higher the average total mass in each point remains 
constant for all t . This seems surprising since a symmetric random walk is recurrent 
in dimensions 1 and 2, but it follows by a rescaling argument and the fact that a 
Brownian motion is point recurrent only in dimension 1 . In dimension 1 we give a 
representation of the asymptotic mass that depends on 

a := K/p 

but not on the strength of the potential y. Let 



m x {t) := (u(t,x)) , (t,x) e [0,o°) x Z d ', 
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denote the expected mass at time t in the lattice point x. The main results of this 
section are for d = 1 , 



Theorem 3. For all x £ Z, 



lim m x (t ) = 1 I ds — — 



i 

1 



as 2 1 



(14*) 



and for higher dimensions 

Theorem 4. For d>2 and all x £ Z d , 

limm T (f) = 1. 

Remark 2. Even though the formula in Theorem [3] looks quite clumsy we find that 
lim t ^ooin x (t) is decreasing in a. It tends to 1/2 as a tends to zero and it tends to zero 
as a tends to infinity. 

The third section is dedicated to analysing the leading order asymptotics of moments 
of the PAM solution from the perspective of the catalyst, i.e., we consider y > and 
the expression u(t,x) := u (t,Y t +x). For p £ N and x = (xi , . . . ,x p ) € V d we denote 
by 

\ i=i 

the p-th mixed moment at x. Moreover, introduce the p-th Hamilton operator on 

r {W d ) by 

J4? p :=£? P + YV P 

where the potential V p is defined as (V p f)(x) = Y%=\ £o(*;)/( x )> an d ^ p acts on 
r (Z pd ) as 

*/Pf(x) = K J £ (f(x + e)-f(x))+p £ (f( Xl+e ,...,x p + e)-f(x)). 

\e\=\ \e\ = \ 

Here, the first term represents the random movement of a collection of p indepen- 
dent random walks accounting for particle diffusion, and the second term arises 
from the shift by the position of the catalyst. By application of the well-established 
Feynman-Kac formula and calculating the generator of the resulting semigroup, we 
obtain the operator representation 

<(/)= (e'^'l) (x), x£Z pd . (2) 

This gives the connection between large time moment asymptotics and spectral anal- 
ysis of the above Hamiltonian. Let us denote by l p the supremum of the / 2 -spectrum 
of Jtf? p . Gartner and Heydenreich |GH06| have shown that, for all p G N and inde- 
pendently of x £ Z d , 
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Hm(l/t)log(u(f J Jc)')=V 

This limit is called p-th Lyapunov exponent. It can be shown by similar methods 
that just as well 

lim(l//)log<(/) =k p , xeZ pcl . 

However, this does not enable us to derive large time asymptotics up to equivalence. 
Assuming the existence of an eigenfunction (v p ) corresponding to X p with certain 
properties, we could on a heuristic level decompose the right hand side of equation 
© as 

ml (t) = e^(l, v p ) ;2 v p (x) +o(e'H x e V d . 
Our next main result contains criteria under which this is indeed possible. 
Theorem 5. Fix K > 0, p > and let one of the following conditions be satisfied: 

(i) p = 1 or p = 2, 7 large enough to ensure A p > 0, 

(ii) p e N, 7 > Ad (Kp + p ). 

Then, there exists a strictly positive and summable I 2 -eigenfunction v p of JtffP cor- 
responding to X p > 0. Assuming v p to be normed in I 2 (Z p ^), the large time asymp- 
totics of the p-th moment are given by 

ml (/)~eVvp(x) llvpHj, f->°o, (3) 
where ||-||j denotes the norm in l l (Z prf ). 

Remark 3. In the case p = l, A p is strictly positive if and only if 1/7 < G K+p (0). In 
this case, the existence of a suitable eigenfunction has been known for quite a while, 
see e.g. MCM94II or IGdH06l . 

Remark 4. For the cases p > 2, the condition l/y< pG K + p (0) is sufficient to have 
positive exponential growth (i.e., A p > 0). The condition 7 > Ad (Kp + p) also im- 
plies exponential growth of the p-th moment. 



2 Moving trap 

This section is devoted to the case 7 < 0. Our main proof tool is the Feynman-Kac 
representation of the solution u given by 

u(t,x) =Ef exp^y 8y,_ s {X s ) djjuo^), (t,x) e [0,-) x Z d . 

Ej, and E? denote the expectation of a random walk with generator kA, pA 
and (fc + p)A, respectively. The subscript z indicates the starting point and the cor- 
responding probability measures will be denoted by P',. By 



6 Adrian Schnitzler and Tilman Wolff 

p t (z) = P£ (X t/K = z ) = Pf (X,/ K = 0) 
we denote the transition probability of a random walk with generator A . 



2.1 Localized initial condition 



In this section we prove Theorems Q] and [2] With the help of the Feynman-Kac 
representation and a time reversal we find that, for all t > and z € Z d , 

M z (t) = Ef El exp \ y f So {X s - Y s ) ds \ = Ef exp \ y [ 5q (Z s ) ds \ . 



2.1.1 Dimensions 1 and 2 



We start with the dimensions where the random walk is recurrent. 

Proof (TheoremU). Using the semi-group representation of the resolvent 
(A - (K + p)Ay l we find that 




This implies, for all A > 0, 

jdte- l 'M (t) = U(l-yjdte- l ' P{K+p)t (0))\ . 
o V o / 

Now the claim for z = follows by a standard Tauberian theorem. The case z 7^ 
follows due to the recurrence of Z. □ 



2.1.2 Dimensions 3 and higher 



A Tauberian theorem is not applicable in transient dimensions because here the 
expected number of particles does not converge to zero. 

Proof (Theorem^}. Let 
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v(z) := HmM-(r) = E? exp \ y f <5o (Z. s ) ds L z g 1 d . 



Notice that the Green's function G K+p is finite in transient dimensions and admits 
the following probabilistic representation. 



Gr+pW =Ef J So(Z s )ds, zeZ d 



That implies v(z) G (0, 1) for all z G Z . Furthermore, we find that v is the unique 
solution to following boundary problem 

UK + p)Av(z) + y8o{z)v(z) = 0, zeZ d , 
1 limui^ v(z) = 1. 



Hence, for all z G Z c ', 



□ 



V(z) = 1 + 7 - -r Gl (Z) 

K+p-yGi(0) 



2.2 Homogeneous initial condition 

In this section we prove Theorems [3] and |4] For a homogeneous initial condition the 
Feynman-Kac representation yields, for all t > and x G Z d , 

m x (t) = £ Ef E^ exp \ y f 5q (X s - Y s ) ds \ So (Y t ) . 
yez" { i J 



2.2.1 Dimension 1 

Let T := inf{/ > 0: X, = Y t } = inf{f > 0: Z, = 0} be the first hitting time ofX and 
Y. The density of T with respect to Pf, z ^ 0, will be denoted by /f. To prove 
Theorem|3] we split m x (t) into two parts 

m x (t):= £ E^E ) y l T> ^(T r ), 
where X and 7 have not met up to time t , and 
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,(»):= £ E ? E v ^<t ex P { 7 / * (X, - y,) ds 1 $, (Y, 



where they have already met by time t. The next proposition shows that m x is asymp- 
totically negligible. Notice that this implies that there is no difference between the 
hard trap (y= — °o) and the soft trap (ye (— °°,0)) case because m x does not depend 
on y. 



Proposition 1. For all x eZ, 



limra T (f) = 0. 



Proof. We can assume without loss of generality that x = since Z is recurrent. Let 
a := inf {t > 0: Z, 7^ Zq} be the the first jumping time of Z. Furthermore, for t > 
let 

*(*):= Eg exp I r| 60 (Z. v ) dj | & (Z r ) . 

In a first step we give an upper bound for the rate of decay of w. Let us abbreviate 
a := 2{k + p). Using the strong Markov property of Z we find 



w(t) = E^l (J> ,e^+Efl a < r e^ 
t 

= e (r-a)t + a f d^e^-^Efl^ 



E| a exp| 7 J 5o(Z H )d M |^(Z,_. v ) 

E^exp{ 7 / do(Z u )du\do(Z t _ s _ r ) 







l-E(t) 



a 



a-y 



(y*f])*w){t). 



Here E denotes the distribution function of an exponentially distributed random 
variable with parameter a — y and y denotes the corresponding density. By iteration 
we find that, for any k > 1 , 



= (!-£)*£ 



a 



a 



a-y 



(*+i) 



Since there exists Ci > such that < Ci (1 + 1) 3//2 for all z ^ and f > 0, we 
see that asymptotically 



< (1 -£(*)) * (d (1 +f)" 3/2 ) = C 2 (1 +r)- 3 / 2 , 
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where C2 is a positive constant. Let :=X — Y andZ^ 2 ^ :=X + Y. Then it follows 
by Holder's inequality that 



«o(t)= I E z ;;,Ef' l T <,exp - Y J 8o(Z^)ds\8 : \zr)8 z (z} 



(2) 



< £ E^'l^exp. 

z,yez V 



2/3 



1/3 



For f > let 



h{t) : E§ 1 e\p 



8 Q (z^) ds \8 Q (z^) 



Obviously h admits the same asymptotic behaviour as w. Fix K > 0. The strong 
Markov property and the central limit theorem yield 



mo(0 < £ 



1 



l?<a-j h(t-s-r) 



2/3 



■(P2(*:+p)«Cy,z)) 



1/3 



Ef ( " W s 



ft (f — i — r) 

/'2( ) c + p) J (y,z)) 1/3 <c^- 1 / 6 '^o 



2/3 



Here C is a positive constant. This proves the claim . □ 

Now we show what m x asymptotically looks like. Recall that a = K/p. 

Proposition 2. For all x G Z, 

1 



1 I' yj(l+a)(l-s)s + g- a 
hmm_ v (f) = 1 / di- 



ll 



as 2 (l 



(l+a) 



Proof. Because of the strong Markov property of X and Y we find 

m x {t) = ^EX5o(i-r)-E E ? E vlT<,^m) 

yez yeZ 

= 1 - £ E*E£l^[E£l£ W-*)]*=t 

yeZ 
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= l-2>XWp(r-T)(r T ). 
It follows by Donsker's invariance principle that 

Mm £ EfEjl T < fjPp(f _ T) (F T ) = / dyE^E^l m p {G ^ {w) (w { \ 



Here and denote two independent Brownian motions that start in the 
origin with variance 2k and 2p, respectively. Their expectations are denoted by 
E}f (1) and E^ (2) , respectively. Moreover, T, W := inf{f > 0: W t W - w/ 2) = y} and 
pf> denotes a Gaussian density with variance 2s. 

Indeed, the application of Donsker's invariance principle is not trivial because we 
have to sum over all where it cannot be applied uniformly. 

LetW"H :=ffW-ff( 2 ),wW :=f( 1 ' + |ff( 2 ' and t| _) :=inf{f >0: W, H =y}. 
Notice that and W i+) are independent. It follows that 

/^ ,,E » ,2, v,^»(^) 

dyEo^'lH p (G) (y) 



' (*+P) 2 



1 

ds ■ 



- s v^> 



^+p)p(l-,),+g 



Now the claim follows by substituting a. □ 
Theorem|3]follows immediately from Propositions[T]and|2] 



2.2.2 Dimensions 2 and higher 



In dimensions 2 and higher, we find that asymptotically the expected mass remains 
constant because a Brownian motion is point recurrent only in dimension 1 . 
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Proof (Theorem^. Let zf' := inf{f > 0: Z, eB £ (Q)} be the first time that the 
process Z hits the centered ball B £ (0) with radius e > 0, and let 

m x (t):= £ E?l£l cz, 8o(Y t ). 

Similarly as in the case d = 1 we find with the help of Donsker's invariance principle 
that 

lim 1 -m () (t) = Km J dxF™ (t^ < . 
However, for d > 2 and x ^ 0, 

limpf ( T f ) < i) =pf ^n{^° ^ 5}) = F " (c < 5) = 0. 

Hence, it follows by monotone convergence that lim,^oomo(f) = 1 which implies 
that lim f _ ) .oo m x (t) = 1 for all xelf 1 . □ 



3 Moving catalyst 

In this section we stick to the homogeneous initial condition mq = 1 and examine 
the case of a randomly moving catalyst, i.e., we consider 7 > 0. 



3.1 Spectral properties of higher-order Anderson Hamiltonians 

Throughout this section, we write k p := super (Jif p ) for all p £ N. Considering the 
first Hamilton operator given by 

J? 1 := (K + p)A + y8 Q , 

the existence of an eigenfunction vi £ I 2 (Z d ) corresponding to its largest spectral 
value, provided that this value is greater than zero, has been widely known for some 
time. The following theorem extends this to the case p = 2 and constitutes the main 
statement of this section: 

Theorem 6. Assume A2 > 0. Then, X2 is isolated in the point spectrum of ffi' 1 
with one-dimensional eigenspace. The corresponding eigenfunction may be chosen 
strictly positive. 

For a start, we restrict the operator to the subspace of component-wise symmetric 
functions 
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© 2 := {/ g I 2 (l? d ) \f(x,y) = f(y,x) Vx,y e Z 2d } , 

which is obviously closed in I 2 (Z 2d ) . Recall the definition of the operators si 2 and 
V 2 from Section[T]and define .e/ 2 and V 2 as their restrictions on the set & 2 above. In 
the same manner, we denote by J^ 2 the restricted second-order Hamilton operator. 
The reader may easily retrace these operators are endomorphisms on & 2 . In partic- 
ular, Jf 2 is a self-adjoint operator on the Hilbert space 6 2 , and it is essential that 
the supremum of its spectrum coincides with A2, which can be shown elementarily. 
Each eigenfunction of jft 2 corresponding to A2 is an eigenfunction of J^ 2 as well. 
Moreover, we expect that an eigenfunction of J$? 2 is, or at least could be chosen 
as, an element of © 2 . In view of that, passing over to & 2 is just a natural approach. 
In the next step, we write g/ 2 + yV 2 rather than J#" 2 in order to emphasize the de- 
pendence on the potential parameter 7, and we establish a further translation of the 
main task: 

Lemma 1. Suppose X > 0. Then, the resolvent operator Rx := (A — g/ 2 )^ 1 exists 
on & 2 , and for all y > 0, we have 

(i) 

Xea(J 2 + yV 2 ) 

(ii) 

X = supa^ + yV 2 ) 
Moreover, for v € & 2 and 7 > 0, 
C Hi) 

{J 2 + yV 2 )v = Xv {R x V 2 )v=y. 

Proof. A Fourier transform reveals that the spectrum of srf 2 is concentrated on the 
negative half-axis, thus (X — ,c/ 2 )~ l exists on I 2 (Z 2d ) for all X > 0. In particular, 
it exists on © 2 , and then it coincides with (A — £/ 2 )~ l as srf 2 is an endomorphism 
on & 2 . Assertions (i) and (iii) follow by rearranging the equations considered and 
applying the resolvent operator. The second relation is shown using the Rayleigh- 
Ritz formula. □ 

As a next step, we introduce an operator 7\ on I 2 CZ d ) having the same spectrum and 
the same point spectrum as RxV 2 and that admits the decomposition 7\ = f W + 
f 2 ). Here, f^ 1 ' is compact and the supremum of ff(f ™) is strictly smaller than the 
supremum of (J (fx). Then, we use Weyl's theorem to obtain that the largest value 
in o(f x ) belongs to the point spectrum a p (fx). The resolvent R x := (A — si 2 ) 
admits the representation 

(Rxf) (xi,xz) = r [ 2 \yi- x uy2-x 2 )f(yi,y2), x u x 2 eZ d , 



y- x ea(RxV 2 ) 1 



y- 1 = super (RxV 2 ) 
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where the resolvent kernel r, : Z 2d 



dfe A 'P (Z f = {xi,x 2 )), x h x 2 G 



Here, Z is a random walk on Z 2rf with generator sd 2 . Then we obtain 

(/? A F 2 /) (jci,JC 2 ) = E r< 2) (y 1 -x 1 ,y 2 -x 2 )[^(yi) + ^(y 2 )]/(3'i,3'2). (4) 

yi,y 2 eZ d 

If we assume / G © 2 , we get 



(^ 2 /)(*i,*a)= E [ r i Cy-'i.-'a) + rjf J (-*i ) y-* 2 )J /Cy,o), 



forxi,JC2 G Z , and in particular 



(/^V 2 /) M) = E [^(y-x,0)+r) L 2, (-x,y)j/(y,0), 
forx G Z rf . Let us therefore introduce the operator 7\ = T^> + on I 2 (Z rf ) with 



f «/» := E ^ 2) (-x,y)/(y), ? (2, /» := E 'Cy-^/ty), * g 2 



Both operators are apparently self-adjoint. The lemma below identifies the spectra 
and point spectra of 7\ and /s^V 2 : 

Lemma 2. For a// A > 0, 

a (f A ) = a (^V 2 ) , o p (f K ) = a p (fl A V 2 ) . 
Proof. The crucial and least intuitive part is to show that 



\l — Tx surjective =>■ jx —R^V 2 surjective. 



(5) 



All other implications are rather straightforward and we omit them for the sake of 
conciseness. Assume jj, — T% is surjective and choose g G © 2 . Define g(x) := g(x,0) 
for xGZ d . There exists / G I 2 (%> d ) with (n — f^)f — gby assumption. We define 



f(x\,x 2 ) := f(xi)So(x 2 ) + f(x 2 )8o(xi) - 8o(xi)8o(x 2 )f(0), x u x 2 G Z cl 
and then, for x\ ,x 2 G Z rf , 
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F(xi,x 2 ) :-- 



f(xi,x 2 ), jci=0orx 2 = 0; 

II- 1 (R x V 2 f{xx t x 2 )+g(x u x 2 )), else. 



We realize that F £ 6 2 and proceed by showing that F is the desired function 
satisfying (n~R X V 2 )F = g. Note that R x V 2 F(x 1 ,x 2 ) = R x V 2 f(x 1 ,x 2 ) for all 
x\ ,x 2 £ If 1 . In the first place, we have 

R x ? 2 F{x u 0) = £ [rf ) (y-x 1 ,0)+rf ) (-x 1 ,y)]/(y,0) 
= f x f{x,) = i if{x l )-g{x l ) 

= flF(x 1 ,0)-g{x 1 ,0), Xl eZ d , (6) 

and by symmetry R x V 2 F(Q 1 x 2 ) = llF(0,x 2 ) — g(0,x 2 ) forx 2 £ Ir. Moreover, 

HF{xi,x 2 )-R x V 2 F(x u x 2 ) 
= R x V 2 f(x x ,x 2 )+g(x u x 2 ) - R x V 2 f(xi ,x 2 ) 
= g{x\,x 2 ), x u x 2 E Z d ,x u x 2 7^0. (7) 

Equations © and (0 yield the desired result (fx — R X V 2 )F = g. Thus, we have 
shown (O. □ 

In the next step, we are able to calculate the supremum of the spectrum of f( 2 \ Its 
value is given in terms of the Laplace resolvent kernel defined by 

r\{x):=J dte- l 'F (X t =x), xeZ d , (8) 
with X a random walk on Z d with generator kA . 
Lemma 3. We have supa(f ( 2 >) = ||f (2) || 2 = ^ +P (0). 
Proof. It will be sufficient to show that 

sup{| M |: M ec7(f( 2 >)}=rf p (0). (9) 

The proof involves a Fourier transform (which we denote by of the operator 
f( 2 \ For / G L 2 ([— n, 7t) d ), the transformed operator reads 

f < 2 >/(Z) = (2*)"* £ e 1 ^ £ rf (y-,,0) / dfce"^/^) 



.vtl 



{2n)- d £ e 1 ^ / dfce-^>V» £ e^Vf (*-y,0) 
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= (jrjr-i/) (/) ( £ e^rf ( z ,0) ) , I G [-%,%) d . 

\zez d ) 

Thus, f ( 2 ) is a multiplication operator and the multiplier 

r(/):= £ e^rf ( z ,0), 

is obviously continuous. Hence, its spectrum is just the closure of the range of that 
multiplier. As each of the two components of a random walk on l? d with generator 
si 1 is just a random walk on Til 1 with generator (k + p)A, we have 

sup |?(/)|=?(0)= £ r?( Z ,0)=rf(0), 

le[-7t,7t) d z£ jd 

and equation (O follows by taking into account that the Fourier transform is an 
isometry. □ 

Lemma 4. Suppose X2 = sup a (j^ 2 ). Then, the operator 7L = r' 1 ' + T ' 2 ' /za.s a 
strictly positive eigenfunction v corresponding to its largest spectral value 1 jy. This 
value is isolated in the spectrum. 

Proof. At first, we realize that T W belongs to the trace class as 
£ (f (1 »5 v ,5,) < f d/e-^<°o, 

and therefore ?W is compact. Then, we explain why sup o^f' 1 ' + T^) > sup<r(f ' 2 ^), 
which together with Weyl's theorem (see e.g. IIRS72I ) yields the existence of an 
eigenfunction. In the end, it remains to show that we may choose this eigenfunction 
strictly positive. 

In order to show that sup <j(T^> + T^) > supcr(f ( 2 )), we recall that 

super (7^) = super (/^V 2 ) = i 

by Lemmas [TJ and 12 and the supremum of a{t^) is equal to ^^(O) by Lemma 
[3] Therefore, it suffices to show that 

i>C'(0). (10) 

Let X\ := super (^f 1 ). In case X\ > 0, it is well-known that 
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compare e.g. Carmona and Molchanov [CM94|. Moreover, as Ai and A2 are the ex- 
ponential growth rates of the first and second moment of u (t,x), Holder's inequality 
yields 

thus a fortiori A] < A2. As A 1— > r^ +p (0) is strictly decreasing (see e.g. equation ([H)), 
and we have shown ( fTob for the case Aj > 0. In case Ai = 0, we have 

\>g k+p {o), 

and we arrive at ( fTOt as G K+p (Q) > r^ +p (0) for all A > 0. Weyl's theorem now 
states that 1/y belongs to the discrete spectrum of f^ 2 since T"> is compact. Con- 
sequently, the value 1 /y is isolated in the point spectrum. Finally, we show that a 
corresponding eigenfunction v may be chosen strictly positive. It suffices to show 
that Tfa is positive in the sense that it maps nonnegative, non-zero functions to pos- 
itive functions. Choose a nonnegative function / arbitrarily and assume f(yi) > 
for some y\ € Z d '. Then, for all i£ Z , 

VW > [^(yi -x,o)+ r ^(-x, yi )] f( yi ) > 0. 

Consequently, we may choose v strictly positive, and the proof is complete. □ 
Let us now prove the main result of this section: 

Proof (Theorem®. Let A2 = super (^ 2 ). The preceding lemma states that there 
exists a strictly positive function v£l 2 (Z rf ) with 7Vv = (1 /y)v. By Lemma|2l there 
exists v>2 € © 2 with R^.V 2 V2 = (1/ 7)^2, as point spectra of both operators coincide. 
Naturally, V2 is also an eigenfunction of Rx 2 V on I 2 (Z 2rf ) . We easily verify that 

R h V»f > 

for all nonnegative, non-zero / G / 2 (Z 2</ ), thus v>2 may be chosen strictly positive. 
Now Lemma Q] yields that V2 is an eigenfunction of J%? 2 and Jj? 2 corresponding to 
A 2 . 

In order to show that its corresponding eigenspace is one-dimensional, let (w;) i6/ 
represent an orthonormal basis of this eigenspace. The w, are principal eigenfunc- 
tions of R^V P that maps nonnegative, non-zero functions to positive functions. 
Hence we may choose all w ( - strictly positive. As two strictly positive functions 
in I 2 (Z ; ' d ) cannot be orthogonal, it follows that |7| = 1, i.e., the eigenspace corre- 
sponding to A2 is one-dimensional. □ 
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We will additionally need that the largest eigenvalue X2 is isolated in the spectrum 
of jfe 2 in order to describe the asymptotic moment behaviour: 

Lemma 5. The value X2 is isolated in g(J^ 2 ). 

Proof. We know by Lemma [4] that y~ l = suptr(7\ 2 ) is an isolated eigenvalue, so 
there exists 8 > with 

'r l -S,r 1 + ~s]n ( j(T l2 ) = {r 1 }. 

Define now 81 small enough to ensure 

||^-e - T h || 2 < -V 2 for a11 e with < £ < 5i . 

It is quickly verified that this is always possible (e.g. by the mean value theorem). 
We can show with a similar argument that supc(7\) depends continuously on X, 
making it possible to find 82 small enough to satisfy 

sup a (f^-e) - SU P ° [ T h) < 5/2 for a11 £ with < £ < 52. 
If we choose now e < 5i A 82, it follows 

T 1 £ super (f h _ £ ) e [y" 1 - 5/2, y" 1 + 5/2 

and by Theorem [8] below, the interval [y -1 — 5/2, y +5/2] contains exactly one 
element of the spectrum of fx 2 + (f^-e ~ T h) = Tfa-e- Therefore, y -1 e p(7\_ £ ) 
and it follows that X2 — £ £ p(^° 2 ) in the usual way by Lemmas Q] and [2] Thus, X2 
is isolated in a (jft 2 ). □ 

Let us in the following present a sufficient condition for the existence of an eigen- 
function with the desired properties that holds for general p G N: 

Theorem 7. Suppose peN and y > 4d (Kp + p). Then, X p = sup a (,^f p ) is posi- 
tive, isolated in the spectrum and belongs to the point spectrum with one-dimensional 
eigenspace. The corresponding eigenfunction may be chosen strictly positive. 

The proof relies on the following theorem from perturbation theory of bounded op- 
erators. It describes the behaviour of an isolated eigenvalue under a bounded pertur- 
bation that is sufficiently small in a certain sense, see Birman and Solomjak MBS 801 . 
Ch. 9.4 for a proof. 

Theorem 8. Let T 1 S denote two self-adjoint operators on a Hilbert space. Suppose 
/J, £ (J p (T) with multiplicity r and 

[ju -e, ju + e]n a (r) = {Ai} 



for some £ > 0. Moreover, assume O (S) C [8\ , 52] for some 8\ < 82 S K with 82 — 
81 < £. Then, the set 
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\ii + Si,n + Sz]r\a{T+S) 

contains only isolated eigenvalues ofT + S whose sum of multiplicities equals r. 

Proof (Theorem^. We have Jtf? p = stf p + yV p , and the idea is to understand the 
generator srf p as a perturbation of the potential yV p . With increasing y, the per- 
turbation sd v remains relatively small, which allows an application of Lemma [8] 
As yV p is a multiplication, its spectrum coincides with the essential range of the 
multiplier and we easily verify that yp is the largest eigenvalue of yV p and has 
one-dimensional eigenspace. Moreover, 

a (yV p ) n [y P -y,yp + y} = { yp} , 

and we may show by a Fourier transform that 

ffK P ) C [-4d(Kp + p),0]. 
Theorem[8]now yields that the set 

e{jr p )n[yp-y,y P ] 

contains exactly one element, which is an eigenvalue with multiplicity one. This 
element must be X p = sup a (Jf? p ) due to the nonpositive definiteness of srf p '. It re- 
mains to show that the corresponding eigenfunction may be chosen strictly positive. 
To that purpose, we consider that v p is also an eigenfunction of e correspond- 
ing to its largest eigenvalue eV Employing the Feynman-Kac representation of this 
operator, we see that it maps nonnegative, non-zero functions to strictly positive 
functions. This means that all principal eigenfunctions are either strictly positive or 
strictly negative. □ 



3.2 Application to annealed higher moment asymptotics 

A natural approach to more exact asymptotics of mixed moments, and the main idea 
proving Theorem[5] is to decompose the semigroup representation 

ag(t) = (tf jeP i)(x). 

Certainly we must consider the initial condition 1 as an appropriate limit of l 2 - 
functions when attempting a rigorous proof. With {E^\^l <G M} the family of spec- 
tral projectors associated with ,¥f p , the spectral theorem for self-adjoint operators 
yields 

m§(t)=e x P t (l,v p )v p (x) + j h Vd£ M (l)(x) (11) 
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for some e > small enough. Here, X p = super (Jt? p ) must be a positive eigen- 
value with multiplicity one that is isolated in a (Jtf p ), and v p is a strictly posi- 
tive and Z 2 -normed eigenfunction corresponding to X p . Beyond that, we need that 
(l,Vp) < °°. If these requirements are met, we may asymptotically neglect the last 
term in (fTTT ). In order to prove Theorem[5] we need the following two auxiliary lem- 
mas that are of pure technical nature and thus given without a proof. The first one 
enables us to approximate the homogeneous initial condition with / 2 -functions. Let 
Q,' 7 := V d n [-t,t] pd for t > 0. 

Lemma 6. For all x e 7L pd , 

(e^"l Q p)(x)~mf(0, t^oo. 

The second auxiliary lemma ensures that the considered principal eigenfunctions 
are summable: 

Lemma 7. Suppose X p = sup o (Jf? p ) > and there exists a corresponding eigen- 
function v p e l 2 (Z pd ). Then, v p e l l (Z pd ). 

Let us now give a concise proof of the main statement: 

Proof (Theorem^. It suffices to show that, for all x € IP , 

e-V,a£(/)— >||v 1 ,|| 1 v 1 ,(jic) (12) 

as t approaches infinity. The spectral representation of t ,J ^ P yields 



for some e > as X„ is isolated in the spectrum. For t large enough, the / 2 -norm of 
the integral is roughly estimated from above by 



< e" 



V, 



< (2' 2 ) 



pd 

2 ^e-' £ , 



which means we may neglect this term and have 



e " Iq? — MNIi v /» 



t -> oo. 



(13) 



The limit ([T2l now follows by equation 113) . Lemma [6] and the triangle inequality. 
This completes the proof. □ 
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